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Abstract 

The standard method used for tagging 6-hadrons in the DELPHI experiment at 
the CERN LEP Collider is discussed in detail. The main ingredient of 6-tagging 
is the impact parameters of tracks, which relies mostly on the vertex detector. 
Additional information, such as the mass of particles associated to a secondary 
vertex, significantly improves the selection efficiency and the background sup- 
pression. The paper describes various discriminating variables used for the 
tagging and the procedure of their combination. In addition, applications of 
6-tagging to some physics analyses, which depend crucially on the performance 
and reliability of 6-tagging, are described briefly. 



(Accepted by Eur. Phys. J. C) 



ii 



J.Abdallah^e, P.Abreu^'*, W.Adam^-*, T.Adye^**, P.Adzic^^, T.Albrecht^^, T.Alderweireld^, R.Alemany-Fernandez^ 
T.AUmendingeris, P.P.AUport^B, S.Almehed^'', U.Amaldi^o, N.Amapane"^, S-Amato^i, E.AnashkinS^, A.Andreazza^s 
S.Andringa^-*, N.Anjos^*, P.Antilogus^S, W-D.Apcll'', Y.Arnoudl^, S.Ask^'^, B.Asman-*", J.E.Augustin^^ 

A. Augustinus*, P.Baillon*^, A.Ballestrero**^, P.Bambade^^ , R.Barbier^*, D.Bardin^*^, G.Barker^^, A.Baroncelli'*° 
M.Bates^*, M.Battaglia^, M.Baubillier^e, K-H.Becks^^, M.Begalli^, A.Behrmann^s, N.Benekos^^, A.Benvenuti^ 
C. Berates, M.Berggren^^, L.Bcrntzon'16 , D.Bcrtrand^, M.Bcsancon-l^ , N.Bcsson-*^, J.Bibby''*'^^ P.Biffi^^, D.Bloch^ 
M.Blom32, M.Bonesini^o, M.Boonekamp'*^, P.S.L.Booth^s, G.Borisov^s, O.Botner^^, B.Bouquet^^, T.J.V.Bowcock^^ 
LBoyko^*, M.Bracko'*^, P.Branchini*", R.Brenner^^, E.Brodet^^, P.Bruckman^o, J.M.Brunet'', L.Bugge^* 
P.BuschmannSe, M.Caccia^^.lO, M.Calvi^O, T.Camporcsi*, V.Canalc^S, F.Carcna*, N.Castro^-*, F.CavalloS, V.Chabaud* 
M.Chapkin-*-*, Ph.CharpentierS, P.Checchia^^, R.Chierici*, P.Chliapnikov'i'i, J.ChudobaS, S.U.ChungS, K.Cieslik^o 
P.Collins*, R.Contrii5, G.Cosme^^, F.Cossutti^o, M.J.Costa^s, F.Couchot^^, B.Crawley^, D.Crennell^*, J.Cuevas^s 

B. D'Almagnc22, J.D'Hondt^, J.Dalmau^"^, T.da Silva", W.Da Silva^^, G. Delia Ricca^O, A.Dc Angclis^O, W.Dc Bocrl^ 

C. DeClercq^, B.DeLotto^", N.De Maria-*^, A.De Min^'', L.de Paula^^, L.Di Ciaccio^^, H.Dijkstra'^, A.Di 3^006"*° 
K.Doroba^^, J.Drees^^'*, M.Dris^^, G.Eigen^, T.Ekelof^^, M.EUert^^, M.Elsing*, M.C.Espirito Santo*, G.Fanourakisi^ 

D. Fassouliotis^-^, M.Fcindt^^, J. Fernandez'*'', A.Ferrer^^, F.Ferro*^, U.Flagmeyer^^, H.Foeth*, E.Fokitis^^ 

F. Fulda-Quenzer22, J.Fuster^a, M.Gandelman^i, CGarcia^s, Ph.Gavillet*, E.Gazis^^, T.Geraligi^, R.Gokieli*'55 
B.Golob''5, J.J.Gomez Cadenas^S'*, G. Gomez- Ceballos^^^ P.Goncalves^*, E-Graziani^o, G.Grosdidier^^, K.Grzelak^s 
J.Guy38, C.Haagl9, A.Hallgren''^^^ K.HamacherSS, K.Hamilton^S, J.Hansen^'*, S. Hangs'!, p.Haulerl^, V.Hedberg^'' 
M.Hennecke^^, J.A. Hernando^^, H.Herr*, J.Heuser^e.*!^ S-CHolmgren-ie, P.J.Holt*, M.A.Houlden^s, K.Hultqvist'"^ 
J.N.Jackson^s, P.Jalocha^", Ch.Jarlskog^'', G.Jarlskog^'', P.Jarry''^, D.Jeans^®, E.K.Johansson*^^ P.D.Johansson*^ 
P.Jonsson^*, C.Joram*, L.Jungermann^®, F.Kapusta^®, M.Karlsson'*^, S.Katsanevas-^*, E.Katsoufis^S, R.Keranen'^^ 

G. Kernel'*^, B.P.Kersevan*''*^, A.Kiiskineni^, B.T.King^s, N.J.Kjaer*, P.Kluit^^, P.Kokkiniasi^, C.Kourkoumelis* 

0. Kouznetsov^*, Z.Krumsteini*, M.Kucharczyk^o, W.Kucewicz^o, J.Kurowska^^, J.Lamsa^, G.Leder^*, F.Ledroit^^ 
L.Lcinonen'"^, R.Leitner^^, J.Lemonne^, V.Lepeltier^^, T.Lesiak^", W-Liebig'"^", D.Likc'^'*, A.Lipniacka'"^, J.H.Lopes^l 
J.M.Lopez^s^ D.Loukas^^, P.Lutz'*^, L.Lyons^**, J.MacNaugliton54, A.Malek^e, S.Maltezos^^, F.Mandl^'*, J.Marco'*^ 
R.Marco*^^ B.Marechal^^, M.Margoni^''', J-C. Marin*, C.Mariotti*, A.Markou^^, C.Martinez-Rivero*^ 
F.Martinez-Vidal^S, J.Masik^'*, N.Mastroyiannopoulosl^, F.Matorras''^, CMatteuzzi^", F.Mazzucato^'^, M.Mazzucato^''' 
R.McNulty25, C.Meroni^s, W.T.Meyer^, E.Migliore'*'^, W.Mitaroff^'*, U.Mjoernmark^'', T.Moa'"', M.Mochi^ 
K.Moenig*'", R.Monge^^, J.Montenegro^^, D.Moraes^^, S.Moreno^", P.Morettinii^, U.Mueller^s, K.Muenich^'^ 
M.Mulders^^, L.Mundim'', W.Murray^*, B.Muryn^^, G.Myatt^"^, T.Myklebust^'l, M.Nassiakoul^, F.Navarria^ 
K.Nawrocki^^, R.Nicolaidou'*^, P.Niezurawski^^, M.Nikolenko^*'^, A.Nomerotski'^'''l'', A. Norman''^, A.Nygren^'' 
A.Oblakowska-Mucha^i, V.Obraztsov**, A.Olshevskii*, A.Onofre^*, R.Orava", K.Osterbergi^, A.Ouraou*^ 
A.Oyanguren^'', M.Paganoni^^, S.Paiano^, J.P.Palacios-^^, H.Palka-^", Th.D.Papadopoulou^S , L.Pape*, C.Parkes-^^ 
F.Parodi^^, U.Parzefall*, A.Passeri"*", O.Passon'"'*', L.Peralta^**, V.Perepelitsa-''-\ A.Perrotta'', A.Petrolini^^, J.Piedra''^ 
L.Pieri*", F.Pierre*^, M.Pimenta^'*, E.Piotto*, T.Podobnik^-^ V.Poireau'*^, M.E.Pol^, G.Polok^O, P.PoropattSO 
V.Pozdniakovl*, N.Pukhaeva^'i^ ^ A.Pullia'^", J.Ramesi-*, L.Ramlerl'^ A.Read'^-*, P.Rcbecclii*, J.Rehnl", D.Reid''*^ 
R.Reinliardt^e, P.Renton^^, F.Richard^^, J.Ridky^"', M.Rivertf's, D.Rodriguez"''^ A.Romero47, P.Roncliese^^ 

E. Rosenbergi, P.Roudeau^^, T.Rovelli^, V.Rulilmann-Kleider''^, D.Ryabtchikov*'', A.Sadovsky^*, L.Salmi^^, J.Salt^^ 
A.Savoy-Navarro^'', U.Schwickerath*, A.Segar^'^, R.Sekulin^*, M.Siebel''"'', A.Sisakian^*, G.Smadja^*, O.Smirnova^'' 

A. Sokolov'*'*, A.Sopczak^s, R.Sosnowski^S, T.Spassov*, M.Stanitzkii'-^, I.Stavitski^^'^s _ A.Stocchi^^, J.Strauss54 

B. Stugu^, M.Szczekowski55, M.Szeptycka^s, T.Szumlak^i, T.Tabarelli*°, A.C.Tafifard^s, F.Tegenfeldt^^ 
J.Timmermans32, N.Tinti^, L.Tkatchevl* , M.Tobin^S, S.Todorovoval'*, A.Tomaradze*, B.Tome^'*, A.Tonazzo^" 
P.Tortosa'^S, P.Travnicek", D.Treille*, W.Trischuk'*^, G.Tristram^, M.TrochimczukSS , CTroncon^^, M-L.Turluer-^^ 

1. A.Tyapkini*, RTyapkin^*, M.Tyndel^*, S.Tzamariasi^, V.Uvarov'''*, G-Valenti^, P.Van Dam^^, J.Van Eldik* 
A. Van Lysebetten^, N.van Remortel^, I. Van Vulpen^-^, G.Vegni-^^, F.Veloso^'', W.Venus^*, F.Verbeuret^^ P.Verdier^* 
V.Verzi^s, D.Vilanova*^ L.Vitale^u, V.Vrba", H.WahlenS^, A.J.Washbrook^s, P.Weilhammer*, CWeiser^^, D.Wicke* 



iii 



J.Wickens^, G.Wilkinson^e, M.Winter", M.Witek^o, O.Yushchenko"''*, A.Zalewska^o, P.Zalewski^s, D.Zavrtanik^s, 
N.I.Ziminis, A.Zintchenkoi®, M.Zupan^^ 

^Department of Physics and Astronomy, Iowa State University, Ames lA 50011-3160, USA 
^Physics Department, Universiteit Antwerpen, Universiteitsplein 1, B-2610 Antwerpen, Belgium 
and IIHE, ULB-VUB, Pleinlaan 2, B-1050 Brussels, Belgium 

and Facultc des Sciences, Univ. de I'Etat Mons, Av. Maistriau 19, B-7000 Mons, Belgium 
■^Physics Laboratory, University of Athens, Solonos Str. 104, GR-10680 Athens, Greece 
^Department of Physics, University of Bergen, AUcgaten 55, NO-5007 Bergen, Norway 
■''Dipartimento di Fisica, Universita di Bologna and INFN, Via Irnerio 46, IT-40126 Bologna, Italy 
® Centre Brasileiro de Pcsquisas Ffsicas, rua Xavier Sigaud 150, BR-22290 Rio de Janeiro, Brazil 
and Depto. de Ffsiea, Pont. Univ. Catolica, CP. 38071 BR-22453 Rio de Janeiro, Brazil 
and Inst, de Fi'sica, Univ. Estadual do Rio de Janeiro, rua Sao Francisco Xavier 524, Rio de Janeiro, Brazil 
^College de France, Lab. de Physique Corpusculaire, IN2P3-CNRS, FR-75231 Paris Cedex 05, Prance 
8CERN, CH-1211 Geneva 23, Switzerland 

^Institut de Recherches Subatomiques, IN2P3 - CNRS/ULP - BP20, FR-67037 Strasbourg Cedex, Prance 
i°Now at Universita dell'Insubria in Como, Dip.to di Scienze CC.FF.MM' via Vallegio 11, 1-22100 Como, Italy 
i^Now at DESY-Zeuthen, Platanenallee 6, D-15735 Zeuthen, Germany 

i^Institute of Nuclear Physics, N.C.S.R. Demokritos, P.O. Box 60228, GR-15310 Athens, Greece 
i^Now at Fermilab (FNAL), Kirk and Pine Streets, P.O. Box 500, Batavia, IL 60510 

i^FZU, Inst, of Phys. of the C.A.S. High Energy Physics Division, Na Slovance 2, CZ-180 40, Praha 8, Czech RepubUc 
^^Dipartimento di Fisica, Universita di Geneva and INFN, Via Dodecanese 33, IT-16146 Genova, Italy 
^''Institut des Sciences Nuclcaircs, IN2P3-CNRS, Universitc de Grenoble 1, FR-38026 Grenoble Cedex, France 
i^Hclsinki Institute of Physics, HIP, P.O. Box 9, FL00014 Helsinki, Finland 

Joint Institute for Nuclear Research, Dubna, Head Post OfHcc, P.O. Box 79, RU-101 000 Moscow, Russian Federation 
^^Institut fiir Experimentelle Kernphysik, Universitat Karlsruhe, Postfach 6980, DE-76128 Karlsruhe, Germany 
^'^Institute of Nuclear Physics, Ul. Kawiory 26a, PL-30055 Krakow, Poland 

Faculty of Physics and Nuclear Techniques, University of Mining and Metallurgy, PL-30055 Krakow, Poland 
22Universite de Paris-Sud, Lab. de I'Accelerateur Lineaire, IN2P3-CNRS, Bat. 200, FR-91405 Orsay Cedex, Prance 
^^School of Physics and Chemistry, University of Landcaster, Lancaster LAI 4YB, UK 
24LIP, 1ST, FCUL - Av. Elias Garcia, 14-1°, PT-1000 Lisboa Codex, Portugal 
■^^Department of Physics, University of Liverpool, P.O. Box 147, Liverpool L69 3BX, UK 

26LPNHE, IN2P3-CNRS, Univ. Paris VI et VII, Tour 33 (RdC), 4 place Jussieu, FR- 75252 Paris Cedex 05, Prance 
^''Department of Physics, University of Lund, Solvegatan 14, SE-223 63 Lund, Sweden 
^^Universite Claude Bernard de Lyon, IPNL, IN2P3-CNRS, FR-69622 Villeurbanne Cedex, Prance 
^^Dipartimento di Fisica, Universita di Milano and INFN-MILANO, Via Celoria 16, IT-20133 Milan, Italy 
^^Dipartimento di Fisica, Univ. di Milano-Bicocca and INFN-MILANO, Piazza della Scienza 2, IT-20126 Milan, Italy 
31IPNP of MFF, Charles Univ., Areal MFF, V Holesovickach 2, CZ-180 00, Praha 8, Czech Republic 
32NIKHEF, Postbus 41882, NL-1009 DB Amsterdam, The Netherlands 

^''National Technical University, Physics Department, Zografou Campus, GR-15773 Athens, Greece 
^^Physics Department, University of Oslo, Blindern, NO-0316 Oslo, Norway 
35Dpto. Fisica, Univ. Oviedo, Avda. Calvo Sotelo s/n, ES-33007 Oviedo, Spain 
^^Department of Physics, University of Oxford, Keble Road, Oxford OXl 3RH, UK 
^'^Dipartimento di Fisica, Universita di Padova and INFN, Via Marzolo 8, IT-35131 Padua, Italy 
3®Ruthcrford Appleton Laboratory, Chilton, Didcot OXll OQX, UK 

^^Dipartimcnto di Fisica, Universita di Roma II and INFN, Tor Vergata, IT-00173 Rome, Italy 
''ODipartimento di Fisica, Universita di Roma III and INFN, Via della Vasca Navale 84, IT-00146 Rome, Italy 
''^Now at Inst, of Physical and Chemical Research RIKEN, 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan 
■^^DAPNIA/Servicc de Physique des Particules, CEA-Saclay, FR-91191 Gif-sur-Yvettc Cedex, France 
''^Instituto dc Fisica dc Cautabria (CSIC-UC), Avda. los Castros s/n, ES-39006 Santander, Spain 
**Inst. for High Energy Physics, Serpukov P.O. Box 35, Protvino, (Moscow Region), Russian Federation 
*^J. Stefan Institute, Jamova 39, SI-1000 Ljubljana, Slovenia and Laboratory for Astroparticle Physics, 

Nova Gorica Polytechnic, Kostanjeviska 16a, SI-5000 Nova Gorica, Slovenia, 

and Department of Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia 
^^Fysikum, Stockholm University, Box 6730, SE-113 85 Stockholm, Sweden 

^'^Dipaxtimento di Fisica Sperimentale, Universita di Torino and INFN, Via P. Giuria 1, IT-10125 Turin, Italy 
**INFN,Sezione di Torino, and Dipartimento di Fisica Teorica, Universita di Torino, Via P. Giuria 1, 
IT-10125 Turin, Italy 

*^Now at Institute of Particule Physics of Canada, Univ. of Toronto, Toronto, Ontario, Canada M5S1A7 
5" Dipartimento di Fisica, Universita di Trieste and INFN, Via A. Valerie 2, IT-34127 Trieste, Italy 

and Istituto di Fisica, Universita di Udine, IT-33100 Udine, Italy 
SiUniv. Federal do Rio de Janeiro, CP. 68528 Cidade Univ., Ilha do Fundao BR-21945-970 Rio de Janeiro, Brazil 
■''•^Department of Radiation Sciences, University of Uppsala, P.O. Box 535, SE-751 21 Uppsala, Sweden 
5=^IFIC, Valencia-CSIC, and D.F.A.M.N., U. de Valencia, Avda. Dr. Moliner 50, ES-46100 Burjassot (Valencia), Spain 



^^Institut fiir Hochenergiephysik, Osterr. Akad. d. Wissensch., Nikolsdorfergasse 18, AT-1050 Vienna, Austria 
^^Inst. Nuclear Studies and University of Warsaw, Ul. Hoza 69, PL-00681 Warsaw, Poland 
sSFachbereich Physik, University of Wuppertal, Postfach 100 127, DE-42097 Wuppertal, Germany 
t deceased 



1 



1 Introduction 

The study of heavy b- and c- quarks is one of the most interesting subjects in exper- 
imental High Energy Physics, directly related to the verification of the Standard Model 
(SM) and the search for its possible violations. Where these may occur is very model- 
dependent, but it may well be that the third generation particles will provide some impor- 
tant clues to new effects. This is a large part of the motivation for studying 6-quarks at 
LEP, where top-quark pair production is kincmatically inaccessible. It is thus important 
to have algorithms for selecting events with 6-quarks while keeping backgrounds small. Ef- 
ficiency and purity or background rejection are important parameters of these techniques. 
Because searches for such deviations from the SM often involve precision measurements, 
it is crucial to have a very well understood and monitored 6-tagging algorithm. 

A further reason for selecting 6-quarks is the search for the Higgs boson. For the SM 
Higgs with mass of relevance for LEP and the Tevatron, the predominant decay mode is 
to bb pairs. Thus tagging 6-jets provides a valuable means of selecting candidates while 
reducing backgrounds to low levels, thereby enabling searches to achieve high sensitivity. 

In this paper the 6-tagging technique developed for the DELPHI experiment at the 
LEP electron-positron collider is described. LEP ran at ccntrc-of-mass energies aroimd 
the Z (91 GeV) over the period 1989 to 1995, and then at higher energies up to 208 GeV, 
before being turned off in 2000. Much of the technique used here would be applicable, 
with suitable modifications, in other experimental situations. 

The lifetimes of 6-hadrons arc around 1.6 ps. This means that flight distances arc of 
order 3 mm for a 35 GeV 6-hadron, this being a typical energy in a 2-jet event at the Z, 
or in a 4-jet event at LEP2 energies. Correspondingly the decay tracks from a 6-hadron 
have non-zero impact parameters^, i.e. when extrapolated backward in space they do not 
pass exactly through the beam interaction region. The scale of these impact parameters 
is cr ^ 400 /im. This is to be compared with the DELPHI experimental resolution a of 
about 

(T = 27©63/(psin^/2^)/im (1) 

where p and 6 are the momentum (in GeV/c) and the polar angle of the track. The 
symbol © denotes the quadratic sum of terms. Eqn. (1) is for the impact parameter 
(IP) in the plane perpendicular to the beam; along the beam direction, the resolution is 
slightly worse. Because the micro-vertex detector is crucial for achieving this accuracy in 
IP measurements, it is described in Section 2. 

The impact parameters provide the main variable for 6-tagging. For all the tracks in a 
jet, the observed impact parameters and resolutions are combined into a single variable, 
the lifetime probability, which measures the consistency with the hypothesis that all 
tracks come directly from the primary vertex. For events without long-lived particles, 
this variable should be uniformly distributed between zero and unity. In contrast, for 
6-jets it has predominantly small values. Details of how this variable is constructed are 
elaborated in Section 3. 

Other features of the event are also sensitive to 6-quarks, and some of them are also 
used together with the IP information to construct a 'combined tag'. For example, b- 
hadrons have a 10% probability of decaying to electrons or muons, and these often have a 
transverse momentum with respect to the 6-jet axis of around 1 GeV/c or larger. On its 
own, the high-pj^ lepton tag would have too low an efficiency for many 6-quark studies, but 
the presence of such a lepton is useful information to combine with the IP measurements. 
The combined tag also makes use of other variables which have significantly different 

^See Section 3.1 for more detailed definitions and discussion of impact parameters. 
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distributions for 6-quark and for other events, e.g. the charged particle rapidities with 
respect to the jet axis. Further details on these variables and the way in which they 
are combined are given in Section 4. The combined tag including the lifetime probability 
and secondary vertex mass, rapidities and fractional energy (described in Section 4.3) was 
used for the measurements at the Z (see Sections 6.1 - 6.4). For most LEP2 6-tagging 
analyses, the transverse momentum missing at the secondary vertex and the transverse 
momentum with respect to the jet axis of any electron or muon were also used in the 
combination. 

The combination method used is optimal for uncorrelated variables. The extent to 
which it is possible to improve on the 'combined tag', for example by using extra infor- 
mation such as the jet energy, is investigated in Section 4.6. The resulting 'equalised tag' 
was used in the Higgs search at LEP2 (see Section 6.5). 

Section 5 contains some technical aspects of the 6-tagging. In particular it describes 
some modifications that were required to the physics generators of the Monte Carlo 
simulation. 

Finally, some physics studies for which 6-tagging plays a crucial role are outlined in 

Section 6. First there is the measurement of the fraction of hadronic Z decays which 
contain 6-quarks (sec Section 6.1). A precise measurement of this quantity requires high 
efficiency tagging, while keeping down the backgrounds from other quarks in order to 
reduce the systematic errors. This is followed by applications of 6-tagging to the mea- 
surement of the production rate of events with 4 6-jets, and of the 6-hadron charged decay 
multiplicity. Section 6.4 describes a measurement of the 6-fraction in 3-jet events, which 
is sensitive to the mass of the 6-quark. This uses anti 6-tagging to select light quark 
events. Finally the crucial reliance on 6-tagging of the search for the Higgs is described 
in Section 6.5. 

2 The DELPHI Vertex Detector 
2.1 Overview 

The silicon vertex detectors of the DELPHI experiment have undergone various up- 
grades throughout the lifetime of the experiment. For the statistics collected in 1991-1993 
it provided measurements in the transverse (i?0) planc^ only [1]. 

The DELPHI Double Sided Vertex Detector (DSVD) [2] was installed in the exper- 
iment in early 1994 and by the end of the Z running at LEP had contributed to the 
reconstruction and analysis of approximately 2 million Z decays. Two of its three layers 
were equipped with double-sided orthogonal readouts, thereby upgrading the IP and ver- 
texing capabilities by adding information from the longitudinal [Rz) plane. The extra 
coordinate helps to associate tracks to vertices where the single i?0 view alone might 
have ambiguities. This upgrade led to about a 30% improvement in the 6-tagging effi- 
ciency at fixed purity. The geometrical layout of the DSVD is shown in figure 1. The 
three layers, termed Closer, Inner, and Outer, were at average radii of 6.3, 9.0 and 10.9 
cm respectively, with the Outer and Closer layers instrumented with the double-sided 
orthogonal readout. The three-layer polar angular coverage was between 44° and 136°, 
with the Closer Layer providing additional coverage in anticipation of the subsequent SiT 
upgrade described below. The transverse view displays the large degree of overlap (up to 

'^DELPHI uses a cylindrical polar co-ordinate system, with the z axis along the beam direction (and the magnetic field 
axis). R and arc the radial and azimuthal co-oidiuatcs in the transverse plane, Q is the polar angle with respect to the 
beam axis. The Cartesian co-ordinates x and y are horizontal and vertical respectively. 
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20% of the sensitive region in the Inner Layer), which was an important ingredient for the 
ahgnment. The average thickness of each sihcon module was 0.5% of a radiation length. 
The z readout was routed via an integrated double metal layer, thus adding negligible 
extra material in the barrel region, and helping to keep multiple scattering to a minimum. 

The DELPHI Silicon Tracker (SiT) [3] was a further upgrade for the physics require- 
ments at LEP2, and the barrel part relevant for 6-tagging was fully installed in 1996. 
Physics objectives of LEP2, such as the measurement of four-fermion processes and the 
searches for the Higgs boson or for super-symmetric particles, required a larger polar 
angle coverage than at LEPl. The design goal was to achieve an equivalent 6-tagging 
performance to the DELPHI DSVD, and in addition to extend this to around 25° in ^, 
after which the 6-tagging capabilities were limited by multiple scattering in the beam- 
pipe. The SiT also incorporated end-caps of mini-strip and pixel detectors for tracking 
in the forward region [3,4]. The geometrical layout of the SiT is shown in figure 2. The 
radii of the layers were similar to the DSVD, but the Outer and Inner layers were doubled 
in length to provide the extra angular coverage. The Closer Layer was double-sided, the 
Inner Layer was double-sided for 21° < 9 < 44° (and the corresponding backward region) 
and single-sided in the centre, and the Outer Layer provided Rcf) and Rz measurements 
from its crossed detector arrangement [3]. The impact of this detector on 6-tagging is 
shown in figure 3. 

2.2 Alignment and Performance 

6-tagging quality relies on excellent alignment of the vertex detector. The starting 
point of the alignment was the information from an optical and mechanical survey before 
installation. This was refined with the information from tracks from Z decays, using a 
stand-alone procedure where the momentum of the track was the only information taken 
from the rest of the DELPHI detector. The precision of the vertex detector hits has 
allowed a number of important effects to be identified, including some common to all 
LEP vertex detectors and certain previously unmeasured properties of silicon detectors. 
They include: 

• coherent deformations, such as a torsion or shear of the entire structure; 

• bowing of the silicon modules due to the different response of the silicon and the 
module support to changes in temperature and humidity; 

• barycentric shift effects, whereby the centres of gravity of the charge clouds of elec- 
trons and holes in the sihcon do not correspond to the mid-plane of the detector, 
nor to each other; 

• acoUinearity of the LEP electron and positron beams leading to Icpton pairs from Z 
decays which cannot be assumed to be back-to-back in the alignment procedure. 

More details can be found in [5]. The precise vertex detector alignment has also led to 
better understanding of other detectors, such as the TPC, the track distortions of which 
were corrected. 

The ultimate performance of the vertex detector with respect to 6-tagging can be 
measured by the IP resolution, which in the Rcf) plane can be parametrised by: 

63 

(yB4 = 27 A*m (2) 

p smafc* 

with p in GeV/c. The IP resolution in the Rz plane for two typical 9 regions can be 
parametrised by: 



4 



Or, = 39 © — /^m (for 80° < ^ < 90°) (3) 
151 

(jRz = 96 ® //m (for 45° < ^ < 55°), (4) 

p 

These equations are the quadratic sums of a constant and a momentum dependent term, 
corresponding to the intrinsic resolution and to the multiple scattering contributions re- 
spectively. For tracks coming from 6-decays, these contributions are of similar magnitude. 
Typical distributions of the IP resolutions as functions of momenta are shown in figure 4. 



3 Lifetime Tagging 

6-hadrons in many aspects are significantly different from all other particles. They 
have a long lifetime, large mass, high decay multiplicity, substantial leptonic branching 
rate, etc. The most important property for the selection of 6-hadrons is their lifetime. 
Among the main features of lifetime tagging are a simple and transparent definition and 
ease of control, since it relies on a single measured quantity, the track IP. 

In this section the definition of the main elements entering in the lifetime tagging 
together with the principles of its construction are given. This tagging itself provides 
efficient separation of the 6-quark from other flavours, which is further enhanced by 
including additional variables (see Section 4). The method of lifetime tagging used by 
DELPHI was originally proposed by the ALEPH collaboration [6]. 

3.1 Impact Parameter 

The general 3-dimensional IP is the minimal distance between the estimated primary 
interaction point and the track trajectory. The decay of a long-lived particle produces 
tracks with large impact parameters, which is not the case for particles from the primary 
interaction. Lifetime tagging is based on this difference. 

For 6-tagging in DELPHI, a slightly different approach is adopted, with a separation 
of the 3-dimensional information into R^) and Rz components. The IP component in 
the i?0 plane is defined as the minimal distance between the primary vertex (PV) and 
the track trajectory projected onto the plane perpendicular to the beam direction. The 
point of the closest approach {Pc) of the track trajectory to the primary vertex in the Rcj) 
plane is also used to define the Rz component of the IP. This is the difference between 
the ^-coordinates of the primary vertex and of the point Pc (see fig. 5). 

According to these definitions, there are two ingredients in the IP computation: the 
parameters of the track trajectory, provided by the track fit, and the position of the 
primary interaction. The parameters of the track trajectory are the track direction given 
by its polar and azimuthal angles (6*, 0) at the point Pq of the closest approach to the 
origin O; and (s^^, Srz)-i the equivalent of the IP components but defined with respect 
to the origin O, rather than with respect to the primary vertex. The reconstruction of 
the primary vertex is explained in the next section. In the approximation that the tracks 
can be regarded as straight lines between Pq and Pc-, the IP components cLr^ and cLrz 

with respect to the primary vertex position V are calculated as: 
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dRct, = Er^- {~e ■ V) 

dRz = + cot e{'u ■V)-Vz 



(5) 



(6) 



Here ll is the unit vector along the track direction in the R(j) plane: ll = {cos (p, sin (p, 0}; 
~e is the unit vector perpendicular to the track direction in R(j) plane: ~e — 

{sin0, — cos(/), 0}; and / = {— cot ^cos 0, — cot ^ sin 0, 1}. Figure 5 illustrates these 
definitions of the IP components. 

The main reason for the separation of the 3-dimensional IP into R(j) and Rz components 
is that the measurement of the particle trajectory in DELPHI is performed independently 
in these two planes with somewhat different precision (see eqns. (2) - (4)). Also the beam- 
spot is smaller in the transverse directions. In addition, there are 3 sensitive layers of 
vertex detector in the R(f) plane and only 2 layers in the Rz plane; the fraction of tracks 
with wrong hit association in the Rz plane is thus higher. The separate treatment of 
the IP components provides the freedom to reject bad measurements in the Rz plane, 
while keeping useful Rcf) information. Finally, the data before 1994 were taken with the 
2-dimensional vertex detector providing track measurements in the R(f) plane only. The 
separate use of the R(f) and Rz information is one of the crucial points of our tagging 
algorithm, significantly infiuencing its structure. 

3.2 Primary Vertex 

The primary vertex is reconstructed for each event using a set of selected tracks and 
the beam-spot position. The beam-spot is the zone of intersection of the two coUiding 
beams of LEP. It has a small size in the Rcf) plane {ax ^ 150/xm, ay less than 10//m), 
while it is several millimetres long along the beam direction. It is relatively stable within 
a fill, and so can be used as a constraint for the primary vertex fit. 

The beam-spot is measured using events which have a vertex formed by at least 3 
tracks with hits in the sihcon strip detectors. These vertices are used to fit the position 
in 3 dimensions and also the x and z size of the interaction region in time periods of 
around 20 minutes. The size of the interaction region in y is not fitted, because it is 
smaller than the corresponding position error, and the value ay = 10 /xm is used. 

The PV position is obtained by minimising the function: 



Here {rf?, ^2} = {d-R^jd^^} is the 2-dimensional vector of IP components for each track 
a entering in the PV fit and Sa is the covariance matrix of the measured quantities 
{^^</>'^Rz}> since measurements in the R(j) and Rz planes are made independently, the 
matrix Sa is almost diagonal. V^^ and a^^ are the beam-spot position and size for the 
X and y coordinates. The first summation in equation (7) runs over all tracks selected 
for the PV fit. Because of our definitions (5-6) of the IP components, the dependence 
of on the vertex position V is quadratic and hence the minimisation of (7) can be 
performed analytically. 

An important part of the PV reconstruction is the selection of tracks and the rejection 
of bad measurements. Tracks with wrong hit associations in the vertex detector, as well 
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as those coming from decays of long-lived particles or from interactions in the detector 
material, bias the fitted PV position and a special rejection procedure attempts to reduce 
this bias. 

For the PV computation, tracks with at least two R(j) measurements and at least 
one Rz measurement are selected. First the fit using all these tracks {Ntr) is performed 
and x^i^Ntr) is computed. After that each track i is consecutively removed and the 
corresponding xf (-/Vtr — 1) is obtained. The track i giving the maximal difference x^i^tr) — 
XiiNtr — 1) is excluded from the fit if this difference exceeds a threshold value A, which 
was set to 6. This procedure is repeated while there are tracks with a difference 
exceeding A. Since the beam-spot position constraint is used for the PV computation, 
all tracks may be rejected for some events. In this case the PV coincides with the beam- 
spot and its covariance matrix corresponds to the beam-spot size. The fraction of such 
events is about 1% for Z hadronic events. 

This fitting procedure gives an average precision of the PV position for qq (where 
q = uds), cc, bb simulated Z hadronic decays of = 36, 44, 60 iim and = 43, 50, 70 /im 
respectively, although the actual precision depends strongly on the number of tracks. 
The somewhat degraded precision for bb events is explained by the smaller multiplicity of 
primary tracks and by tracks from 6-hadron decay occasionally included in the primary 
vertex. 



3.3 Error and Sign of Impact Parameter 

Since the PV position is used in the definition of an IP, the impact parameters of all 
tracks included in the PV fit are correlated with each other; their correlation coefficient 
is about 0.2. From equations (5-7) and the standard error propagation formalism, the 
error on the R(f) IP is given by: 

2 _ / - {^R4,f if the track is included in the PV fit 



" I + otherwise 

with similar equations for a%^. Here cr|^^ {<^rz) is the error on er^ {^Rz) coming from the 
track fit and cr^^ is the error from the PV fit, and includes implicitly the influence of all 
other impact parameters. More explicitly: 

«^)' = J2'^Sr,e, (9) 

i<ir = Y^^rsr^h (10) 

where Sj^ is the covariance matrix of the primary vertex fit, ~e and / are defined in 
Section 3.1, and repeated indices imply summation. The simplicity of the final equations 
is a consequence of our choice of IP components. 

Using equations (5, 6 and 8) the track significances Sr^ and Srz are defined simply 

as: 

Sr,I, — dR^/aRrj) (11) 
Srz — dRz/aRz (12) 



The track significance thus compares the measured value of the IP with its expected 
precision. This quantity is used as an input variable for the lifetime tagging. Tracks from 
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decays of long-lived particles (r's, b-, c- and s-hadrons) often have large IPs, significantly 
exceeding aR^p and ajiz- 

Equations (5-6) define the magnitude of IP components and their geometrical sign, 
while in the 6-tagging method and throughout this paper the lifetime sign for IP is used. 
It requires knowledge of the flight path of the long-lived particle. In the simplest case 
the flight path is approximated by the direction of the jet^ to which the given particle 
belongs. Often the decay point of the long-hved particle can be reconstructed (see Section 
4.2); in this case the flight direction is defined as the direction from the primary to the 
secondary vertex. As can be seen in fig. 6, this improves the measurement of the flight 
direction. The azimuthal angle precision becomes slightly better than that for the polar 
angle because the vertex detector is more precise in the Rcj) plane. To obtain the lifetime 
sign of the Rcf) and Rz IPs, the point of closest approach in space of the track to the 
estimated 5-flight path is computed and the sign is set negative (positive) if this point 
is upstream (downstream) of the PV position. The significance is assigned the same sign 
as the IP. 

With this definition, tracks from decays of long-lived particles have predominantly 
positive signs while tracks coming directly from the PV are equally likely to be positive 
or negative. For 6-tagging, tracks with positive IP are used, thus reducing by half the 
number of background tracks. 

The distributions of positive and negative Re/) significance are shown in fig. 7. The 
excess of positively signed tracks with large significance is clearly seen. 



The distribution of the negative track significance is determined mainly by tracks 
coming from the PV, including scatters in the detector material, tracks with wrong hit 
association etc, while the contribution of tracks coming from decays of long-lived particles 
is about 1%. This distribution can thus be used to define the probability P{S^) for a 
track from the PV to have the measured value of the modulus of its significance exceeding 
the value S'^. This function is obtained by integration of the probability density function 
of the negative significance f{S) from to infinity and assuming that P{S^) is the same 
for primary tracks with either positive or negative significance: 



By definition, tracks from the PV should have a flat distribution of P{S^) between 
and 1, while tracks from decays of long-lived particles and with large positive values of S'^ 
have small values of P{S^), reflecting the small probability for tracks from the primary 
vertex to have such large values of the IP and hence of S^. As an example, fig. 8 shows 
the distribution of P^S*^^) for tracks with positive IP. The peak at small values of PiS*^) 
is produced mainly by the long-lived particles. The transformation from significance to 
track probability is referred to as the calibration of the detector resolution. 

For LEPl analyses, the above calibration was performed using tracks from Z decays. 
At LEP2, there was the possibihty of again using Z data for cahbration; each year, short 
runs at the Z were taken before the start of and also interspersed with the high energy 
running. Alternatively, the calibration could be carried out using the same type of data 
as used to perform the relevant physics analysis. Thus for the Higgs search of Section 
6.5, 4-jet events were also used for cahbration in the channel where the H and Z both 

"^The default jet clustering algorithm is JADE, with j/cut set at 0.01. However, the user of the b-tagging package has 
the option of using any jet algorithm. 



3.4 Track Probability 




(13) 
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decay to 2 jets, while in the corresponding channel where the Z decays to two neutrinos, 
calibration was performed using '2-jets + missing energy' events. The use of calibration 
samples closely related to the data sample in principle allows for the following possible 
effects: 

• because of possible movements of the relative positions of the different parts of the 
vertex detector with respect to the rest of DELPHI, the calibration could be time- 
dependent; 

• track confusion and lifetime-signing (and hence calibration) could depend on the 
event topology; 

• the IP resolution changes with polar angle Q. The jet distribution in Q depends on 

the particular physical process considered; 

• the IP resolution is also energy dependent. The jet and track energy spectra depend 
on the physical process. 

The calibration was performed separately for categories of tracks with different life- 
time sensitivities. The categories were determined by the number of associated VD hits. 
A small number of VD hits associated with a track is often caused by incorrect recon- 
struction, and the significance distribution of such tracks has a larger non-Gaussian tail. 
By using a different track probability for them, this difference was taken into account. 
This approach was also used for the analysis of the data collected in 2000, when one 
out of the 12 sectors of the TPC was not operational during the last part of the data 
taking; tracks reconstructed without the TPC have worse precision, which requires using 
a separate track probability for them. 

Another property of the track probability is that it can be defined directly from the 
data. This is very important, in that it allows the calibration of the detector resolution 
independently of the simulation. Such calibration allows to take into account possible 
differences between data and simulation. As a consequence, it also reduces the systematics 
due to detector effects in physics measurements. 

For the construction of P(5'°) it is important to reduce the contribution of tracks 
coming from the decay of long-lived particles. Using the negative significance distribu- 
tion partially solves this problem. Additional suppression of the lifetime information 
is achieved by applying anti-6 tagging to the event sample used for calibration. This 
anti-6 tagging is based on tracks with positive IPs, and hence does not bias the nega- 
tive significance distribution. The anti-6 tagging reduces the fraction of 66 events in the 
selected sample of hadronic Z events from 21.6% to less than 5% and the contribution 
of tracks from decays of 6-hadrons is reduced correspondingly. Additional selection cri- 
teria for tracks used for calibration decrease the contribution from the decay products 
of light long-lived hadrons [K^ ^ hyperons) and hence reduce the tail of the significance 
distribution. They are specified in the next section. 



3.5 Lifetime Probability 

Track probabilities are directly used to construct a hfetime probabihty [6]. For any 
group of N tracks it is defined as: 

Nr^+Nhz-^ Nr^ Nr^ 

Piv = n • (- n)Vj!, where H = J] P{S'r^) ■ ^(-^k) (14) 

j=0 i=l 1=1 

Here P{S}^^), P{S'^f^^) are the track probabilities and Nr^, Nrz are the number of R(j) 
and Rz IPs used in the tagging. The definition of Pn thus ignores the small off-diagonal 
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elements of the IP error matrix and the correlation between different IPs coming from 
the use of the common PV. 

The variable Pn has a simple and straightforward definition and can be computed 
for any group of tracks (e.g. a jet, hemisphere or whole event) which makes it flexi- 
ble and easily adjustable to different physics applications. No other combination of IP 
measurements was found to give a better selection of 6-quarks. 

An attractive feature of lifetime tagging is that it is constructed using only the track 
IPs. This provides the possibility of achieving a good description of the 6-tagging effi- 
ciency by the accurate tuning of the track resolution in simulation, as described in Section 
3.6. It allows a significant decrease in the systematic uncertainties due to detector effects 
in physics measurements. 

The meaning of the variable Pn is very similar to that of track probabihty P{S): 
it is the probability for tracks coming from the PV to have the product of their 
track probabilities exceeding the observed value. It varies between and 1 and has a flat 
distribution for any group of uncorrelated tracks coming from the PV. The contribution 
of tracks from secondary decays shifts P/v to lower values, producing a peak near 0. 

The flat distribution of Pjv for primary tracks can be verified by computing P^ for 
the sample of all tracks with negative impact parameters in anti-6 tagged events. As 
explained in Section 3.4, the contribution of tracks from decays of long-lived particles in 
such a sample is small. The distribution of is shown as the dotted curve in fig. 9. It is 
relatively flat, although there is a small peak near zero. This peak is produced by tracks 
from decays of long-lived particles, which are occasionally assigned to have negative IPs 
because of the error in their flight direction estimate. However, the vahic of this excess 
is signiflcantly less than the peak in the distribution of P^, computed using positive IPs. 
This distribution is shown as the solid curve in flg. 9. The latter peak is mainly produced 
by the hb events, as can be seen from flg. 9. 

The separation of tracks into two samples depending on the sign of their IP is very 
important. The sample of negative IP tracks is used for the calibration of the detector 
resolution and the quality of this calibration is verifled by the P^ distribution. In contrast, 
positive IP tracks are used for 6-tagging. Thus the samples of tracks used in the calibration 
and in the analysis do not overlap. The P^ distribution also gives a good estimate of the 
background level from light quarks at the corresponding value of P^. 

The P0 and Rz components enter in eqn. (14) separately. As explained in Section 3.1, 
the fraction of wrong measurements in Rz is higher. Therefore, tighter selection criteria 
are applied to tracks for the Rz IP, and for some tracks only the R4> measurement is 
used. Thus, the separate treatment of IP components also allows the use of information 
which would otherwise be lost. 

More speciflcally, the conditions applied to the tracks are as follows. All tracks with 
positive IP and at least one measurement in the VD are candidates for lifetime tagging. 
Tracks coming from reconstructed or A decays^ are rejected. Both P0 and Rz IP 
components are required to be less than 0.2 cm, although this condition is removed if the 
track comes from a reconstructed secondary vertex (see Section 4.2). 

One more parameter is used to provide additional suppression of bad Rz measure- 
ments. It is the distance D of closest approach in 3-dimcnsions between the track and 
the expected flight path of the long-lived particle, deflned in the section 3.3. All tracks, 
both from the primary and secondary vertices, should have a small value of D provided 
the secondary vertex is close to the estimated flight direction. Therefore a large value of 
D with respect to its expected precision ctd is used to identify wrong IP measurements. 

*For reconstruction procedure see [7]. 
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Rz IP measurements are rejected ii D/ao exceeds 2.5. Both Rcj) and Rz measurements 
are excluded from lifetime tagging if D/(Td exceeds 10. 

Figure 10 shows the performance of the lifetime probability applied to simulated 
hadronic decays of the Z. The figure shows the efficiency of the 6-quark selection versus 
the contamination of the selected sample by {u, d, s, c) flavours {Nudsc/ {Nudsc + Nb)). The 
suppression of background flavours is shown for tagging of one jet (i.e. using only tracks 
from the given jet), and for the whole event. Event tagging is more efficient because 
?>-quarks are produced in pairs. 

6-tagging using only lifetime information is rather efficient and is sufficient for the 
needs of many physics applications. An important feature is its simple control of the 
tagging efficiency in simulation. Its performance is however substantially enhanced by 
including additional discriminating variables, such as the mass at the secondary vertex 
or the presence of energetic leptons. This method of 'combined 6-tagging' is described in 
Section 4. 

3.6 Tuning 

Almost all precision measurements and searches for rare processes rely on a comparison 
of the observed data distributions with those predicted by a detailed simulation. For this 
comparison both the generation of the intrinsic physical processes and the simulation of 
detector response must be as realistic as possible. For the selection of events containing 
6-hadrons, the most important variables are the track IPs, therefore the description of the 
IP resolution can significantly influence the physics result and the value of the systematic 
uncertainty. 

The generated events in the DELPHI experiment are processed by the detector sim- 
ulation package [8] and the same reconstruction program [9] as for the data. For the 
simulation, the reconstruction program first applies some additional smearing to the re- 
construction inputs to improve agreement with the particular data set being represented. 
For the vertex detector, this includes applying corrections for inefficient regions, adding 
noise hits, and randomly modifying the positions of the modules to simulate the effects 
of residual misalignments in the real data. 

However, even after this procedure some disagreement between data and simulation in 
the track resolution description remains. This difference can be clearly seen, for example, 
in the distribution of the track significance (see fig. 11). Any disagreement in this 
quantity can result in a large discrepancy in the 6-tagging description. 

A detailed description of the method including the correction of the detector resolution 
in the i?0 plane for the initial micro- vertex detector [1] is given in [10]. The application 
of this method for the tuning in the Rz plane for the DSVD at LEPl and the SiT at 
LEP2 is similar and consists of the following steps: 

• the appropriate parameterisations of the negative lifetime-signed IP {du^ and duz) 
distributions are determined; 

• the numerical coefficients for these parameterisations are extracted from the data; 

• the errors of dR^ and dRz given by the track flt are corrected both in data and in the 
simulation according to the parametrisation obtained while the correlation between 
dR^ and dR^ is not changed; 

• additional smearing of R(j) and Rz IPs in simulation is performed in order to repro- 
duce the observed real data distributions. 

The improvements in the signiflcance description after applying this method can be 
seen in flg. 12. Fig. 13 shows the data to simulation ratio of the selection efficiency as 



11 



a function of the cut on the 6-tagging variables. For the non-tuned version of 6-tagging 
(dashed hne) the difference between data and simulation is very significant for strong 
6-tagging cuts, corresponding to purer samples of B events. The tuning results in better 
agreement for both the lifetime and the combined 6-tagging variables, the latter being 
described in the next section. The remaining differences between data and simulation 
can be explained by the uncertainties of the modelling of B decay and to a lesser extent 
its hadronic production. 

This tuning procedure is incorporated in the 6-tagging package and is used in all 
DELPHI measurements involving 6-quark selection. 



4 Combined Tagging 

Efficient utilisation of different properties of 6-hadrons requires the development of 
a technique for their combination into a single tagging variable. The simplest solution 
of applying some system of cuts on different discriminating variables, which was tried 
in other collaborations [11, 12], is not optimal due to a significant overlap between the 
signal and background for some of them. Instead, DELPHI uses a likelihood ratio method 
of variable combination [13,14]. This approach has the important advantage of being 
technically very simple while at the same time providing powerful separation of signal 
and background. For independent variables, it gives optimal tagging, i.e. the best possible 
background suppression for a given signal efficiency [15]. It can easily be extended to 
any number of discriminating variables, and can deal with different numbers of variables 
in different events. However, its practical application requires the careful selection of 
variables with reduced correlations among them. The description of this likelihood ratio 
method, the set of variables used and the performance of the DELPHI combined 6-tagging 
is given below. 



4.1 Description of Method 

The combined tagging variable y in the likelihood ratio method is defined as: 

_ P^'^jxi, ...,Xn) , . 

^^^^ 

where f^^^{xi, are the probability density functions of the discrim- 

inating variables Xi, for the background and the signal respectively. The selection 

of all events with y < yo gives the optimal tagging of the signal. It should be stressed 
that such tagging is absolutely the best for a given set Xi, ...,Xn of variables. 

In practical applications the determination and utilisation of multi-dimensional prob- 
ability density functions is quite difficult for n > 2. The solution consists in a special 
selection of discriminating variables having reduced correlations among them. In the 
limit of independent variables^, expression (15) becomes: 

y = ny4=n!"; (16) 

i=l Ji \•^^) 1=1 

fi'\x,)/fr{x,) (17) 

^Two variables are independent if, for the signal and for each separately treated background component (e.g. c and 
uds), the distribution of one is independent of any selection on the other. 
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where f-^ (xj), fi^^{xi) are probability density functions of each individual variable Xi for 
the background and signal, and are determined from simulation. 

This scheme is used in DELPHI to construct the combined 6-tagging. For each indi- 
vidual variable Xi the value yi is computed from (17); the combined tag y is defined as the 
product of the yi. It is not exactly optimal any more, because the discriminating variables 
are not independent, but the variables are chosen such that the correlations between them 
are small enough that the resulting tagging is very close to optimal. Furthermore, the 
efficiencies and mis-tag rates are determined from simulation (and sometimes from the 
actual data), thereby taking into account any small correlations. 

In DELPHI all discriminating variables and the 6-tagging itself are computed indepen- 
dently for each jet in an event, where ideally all tracks coming from the fragmentation of 
the 6-quark and from the decay of the 6-hadron are combined in one jet by a jet clustering 
algorithm. In this case the background for the h quark selection can be separated into 
two different parts - jets generated by c-quarks and by light (g = -u, d, s) quarks. These 
two components are independent and have very different distributions of discriminating 
variables. 

To define the extra discriminating variables for the 6-tagging, tracks are selected so 
as to come preferentially from 6-hadron decay. For this purpose all jets in an event are 
classified into 3 categories. In the first category all jets with one or more reconstructed 
secondary vertices are included. A reconstructed secondary vertex provides a clean se- 
lection of 6-hadron decay products and a large number of discriminating variables can 
be defined in this case. If the secondary vertex is not reconstructed, tracks from the B 
decay are selected by requiring the track significance probability to be less than 0.05, and 
the second category includes all jets with at least 2 such offsets. This criterion is less 
strong, allowing more background jets to pass the cut. Finally, if the number of offsets 
is less than 2, the jet is included in the third category and in this case only a reduced 
set of inclusive discriminating variables, like the lifetime probability (see Section 3.5), 
is used. In Z hadronic events the fractions of jets classified into categories 1, 2, 3 are 
44%, 14%, 42% respectively for 6-quark, 8%, 8%, 84% for c-quark and 0.6%, 2.8%, 96.6% 
for light quark jets. 

The tagging variable for a jet of category a is defined as: 

ya = </<YlyU + nMllyl^; (18) 

i i 

where fi^ai^i)^ fi,ai^i)j fi,ai^i) probability density functions of Xi in jet category 

a generated by uds, c and b quarks respectively and n^, and are their normalised 
rates, such that ^^^-^ = Rq, Xl^o = ^c, and Xl^a = ^b- Rq- Rc and Ri, are the 
normalised production rates of different fiavours and Rq + Rc + Rb = ^■ 

As can be seen from eqn. (18), the classification into different categories effectively 
works as an additional discriminating variable with the discrete probabilities given by 
For example, the 6-purity of a sample of jets with reconstructed secondary vertices 
is about 85%. However, the primary purpose of this separation is to allow the use of a 
larger number of discriminating variables when a secondary vertex is found. The search 
for the secondary B decay vertex is thus an important ingredient of DELPHI 6-tagging. 

It is often convenient to define Xj^^t — — logio Va as the jet tagging variable, and this 
variable is used in all applications described in Section 6. The event variable X^v is 
defined as the sum of the largest two Xj^t values for the individual jets in the event. 
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4.2 Secondary Vertex Reconstruction 

A secondary vertex (SV) is searched for in each jet of the event. In the first stage 
all possible combinations of pairs of tracks are selected as SV candidates if they have a 
common vertex with the of the fit less than 4. After that all tracks from the same 
jet are tested one by one for inclusion in a given SV candidate. The track producing 
the smallest change A of the vertex fit is included in the SV candidate if this change 
does not exceed the threshold A = 5. This value and all other numerical parameters 
of the algorithm were selected by optimising the efficiency of the SV reconstruction and 
background suppression. This procedure is repeated until all tracks satisfying the above 
condition are included in the SV candidate. The SV candidate is rejected if the distance 
to the primary vertex divided by its error is less than 4. Additionally, at least two tracks 
in the SV candidate are required to have VD measurements in both R4> and Rz planes. 

The decay of the 6-hadron is usually followed by decays of one or two D mesons, thereby 
producing several secondary vertices. It thus often happens that some secondary tracks 
cannot be fitted to a single secondary vertex. However, the spatial distance between 
any secondary track and the flight trajectory of the 6-hadron should be small since the 
D-mesons tend to travel in the direction of the initial 6-hadron. Using this property 
some tracks of far-decaying D mesons can be recuperated, which is important for the 
computation of such quantities as the 6-hadron mass. The flight trajectory of the b- 
hadron is deflned as the vector from the primary to the secondary vertex. Any track 
from the same jet having an R(j) or Rz component track probability less than 0.03 and 
not included in the SV flt is attached to the SV candidate if its distance in space to the 
flight trajectory divided by its error is less than 3. Although not included in the SV flt, 
such tracks are used in the computation of all discriminating variables. 

Sometimes the decay vertex of a D-mcson can be reconstructed separately from the 
6-hadron decay vertex, which then results in two or more secondary vertices in the same 
jet. In that case, all tracks included in these vertices are combined for the computation 
of the 6-tagging discriminating variables. 

Three additional criteria are used to suppress the background of light quarks in the 
sample of jets with secondary vertices. The flrst one makes use of the momentum vector 
of the 6-hadron. This is deflned as the sum of the momenta of all tracks included in 
the SV candidate. Additionally, the momenta of all other neutral and charged particles 
with pseudo-rapidity exceeding 2 are also included; the pseudo-rapidity is computed with 
respect to the flight direction of the 6-hadron. Then, the trajectory directed along the B- 
momentum and passing through the SV position is constructed and the impact parameter 
Ssv of this trajectory with respect to the primary vertex is computed. For a real 6-hadron 
the momentum direction should be close to the flight direction and Ssv should be small 
compared to its error^ '^^sv while for a false secondary vertex the flight and momentum 
directions are much less correlated. Therefore SV candidates with (Ssv/'^Ssv)'^ > are 
rejected. For the second criterion the lifetime probability using all tracks included in the 
SV candidate is computed and the candidate is rejected if this probability exceeds 0.01. 
The third criterion requires the distance between the primary and secondary vertex to 
be less than 2.5 cm, because the contribution of false secondary vertices and of strange 
particle decays becomes rather high at large distances. Any background jet with a very 
distant SV would give an extremely strong 6-tagging value and this cut effectively rejects 
such cases. 

''crggy is computed using the uncertainties in the positions of primary and secondary vertices and the uncertainty in 
the B-hadron flight direction estimate, obtained from the simulation. 
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Candidates remaining after these selections are considered as reconstructed secondary 
vertices. With this procedure a SV is reconstructed for about 44% of jets with 6-hadrons 
(50% for jets inside the VD acceptance). The &-purity of the sample of jets with a recon- 
structed SV is about 85% for hadronic decays of the Z, which should be compared with 
the initial 6-purity of about 22%. More than one SV in a single jet is allowed, reflecting 
the possibility of cascade {B — > D) decays giving rise to distinguishable secondary ver- 
tices. In this case the tracks from all secondary vertices are combined for the computation 
of the SV discriminating variables. 

4.3 Discriminating Variables 

In this section the discriminating variables used in 6-tagging are described. All defini- 
tions arc given first for jets with reconstructed secondary vertices. Then the modifications 
for other jet categories are described. All discriminating variables, except the transverse 
momentum of a lepton, are computed using the charged particles included in the sec- 
ondary vertex. 

The jet lifetime probability, Pj~, is constructed using equation (14) from the positive 
IPs of all tracks included in the jet. 

The mass of particles^ combined at the secondary vertex, Mg, is very sensitive to the 
quark flavour. The mass at the secondary vertex in a jet generated by a c-quark is limited 
by the mass of the £)-meson, which is about 1.8 GeV/c^, while the mass in a 6-jet can go 
up to 5 GeV/c^. The limit of 1.8 GeV/c^ for the c-jets can be clearly seen in flg. 14(a). 
Some c-jets do have a higher value of Mg due to tracks incorrectly attached to the SV. 

The fraction of the charged jet energy included in the secondary vertex, Xf^, reflects 
the differences in the fragmentation properties of different flavours. The fragmentation 
function for the c-quark is softer than for the 6-quark, as seen in the distribution of Xf^ 
in fig. 14(b). 

The transverse momentum at the secondary vertex. Pi, first introduced by the SLD 
collaboration [12], takes into account missing particles not included in the SV definition. 
P* is defined as the resultant transverse momentum (with respect to the 6-hadron's es- 
timated fiight direction) of all charged particles attached to the SV. Missing particles 
can be neutrinos from scmilcptonic decay, other neutral particles or non-reconstructed 
charged particles. In all cases, due to the high mass of the 6-hadron, the value of P* for 
6-quark jets is higher, as can be seen from fig. 14(c). 

The rapidity of each track included in the secondary vertex, Rf, is quite a strong dis- 
criminating variable, significantly improving the 6-quark selection. Although a 6-hadron 
on average is produced with a higher energy, the rapidities of particles from P-decays 
are less than those from D-meson decay, as can be seen from fig. 14(d). This is mainly 
explained by the higher 6-hadron mass. The variable P*^ is defined for each track in the 
SV and the corresponding variable for each track is used in (18) for the computation of 
the 6-tag. Although there is overlap between the signal and background for an individual 
track rapidity, because of the large number of secondary tracks the inclusion of all the 
rapidities in the 6-tag results in a significant gain. 

The transverse momentum of an identified energetic lepton, p\. It is independent of 
the track IP and can be defined for any category of jet containing a muon or electron. A 
more detailed description of this variable and of some specific features of its inclusion in 
the combined 6-tagging scheme are given in the next section. 

^For computation of discriminating variables, such as mass or track rapidity, elsewhere in this paper, charged particles 
are given the pion mass, whereas neutrals (except for and A) are assumed to be massless. 
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All these variables are defined for the first category of jets (with reconstructed SV). 
For the two other categories, a reduced set of variables is used. For jets with at least 
two offsets, the jet hfetime probability, the effective mass of all tracks with offsets, their 
rapidities and any lepton transverse momentum are computed. For jets with less than 
two offsets the effective mass is not used, as there is no reliable criterion for selecting the 
particles from i?-decay; however the track rapidities for all tracks with positive IP are 
still included in the tagging. The ratios of probability density functions are computed 
separately for each jet category. The possibility of treating in the same way different 
categories of events with different sets of discriminating variables is a very important 
feature of the likelihood ratio method of 6-tagging. 

The distributions of M^, Xf^, P* and Pf are shown in figure 14. These distributions 
are shown for 6-quark jets and also for c-quark jets, the latter constituting the main 
background for 6-tagging. 

Combined 6-tagging using the complete set of discriminating variables performs much 
better than the simple lifetime tagging. This is illustrated in fig. 15. The performance 
is tested using jets of Z hadronic decays. The figure shows the contamination of the 
selected sample by other flavours {Nudsc/ {Nudsc + Nhj) versus the efficiency of 6-jet selec- 
tion. Compared with the tagging using only Pj^ , combined tagging provides much better 
suppression of background, especially in the region of high purity. A very pure sample 
with contamination below 0.5% can be obtained for a sizable b efficiency, which opens 
new possibilities for measurements with 6-hadrons. 

4.4 Lepton Tagging 

Leptons with high transverse momentum have long been used in a variety of ways to 
identify the quark fiavour of the jet from which they originate. This section describes the 
inclusion of this information within the standard DELPHI 6-tagging algorithm. 

Knowing the probabilities pi''^^^ of finding a lepton in a light quark-jet, a c-jet or a 
6-jet respectively, and the transverse momentum distributions f'^''^'^{p\) of these leptons, 
the contribution of identified leptons to the global discriminating variable (18) is: 



The quantities appearing in this expression are extracted from a sample of simulated 
hadronic Z decays, where all reconstructed particles are clustered into jets. 

Reconstructed particles are identified as leptons if they satisfy a tight electron tag 
(from energy loss by ionisation in the TPC, or from associated energy deposits in the 
electro-magnetic calorimeters), a tight muon tag (from the muon chambers only), or a 
standard muon tag confirmed by a minimum ionisation energy deposit in the hadron 
calorimeters. Detailed descriptions of these different categories of tags and of their per- 
formance are given in ref. [7]. The quantities pi^'^'^ are simply defined as the fractions 
of jets of the corresponding flavour inside which a lepton is identifled, and are 3%, 9.8% 
and 18.7% for light quarks, c-quarks and 6-quarks respectively. 

The transverse momentum of the lepton is evaluated with respect to the jet to which 
it belongs, when the lepton momentum is subtracted from the total jet momentum. The 
tagging contribution of the leptons is then obtained from the ratio of the transverse 

momentum distributions /^'^'^•'(p;) and P{p\)- 




(19) 
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Figure 16 shows the superimposed distributions of transverse momenta for leptons 
found in 6-jets, c-jets, and Ught quark-jets. The agreement of these with the data ob- 
tained from real hadronic Z decays (recorded in 1994) is excellent. Because of the small 
branching ratio of 6-hadrons to high leptons, the intrinsic discriminating power of lep- 
ton tagging is weaker than that of lifetime and secondary vertex information. However, 
it does provide useful information since it is fully independent of the above, and supple- 
ments the 6-tagging in particular when no other hints of B decays are found (for example 
in case of fast decay, or decay outside the vertex-detector acceptance). 

In principle, the distribution in transverse momentum at the secondary vertex P* 
depends on whether or not there is a lepton present, because it is accompanied by an 
unseen neutrino. However, because semi-leptonic decays are multi-body, there is not 
a strong correlation between the lepton and the neutrino transverse momenta, and the 
difference in the P* spectra (i.e. with and without charged leptons) is small, as can be 
seen from fig. 17. 

The use of the algorithm at LEP2 requires some additional care. In this case, data 
available for detector calibration are less abundant, and some experimental aspects are 
less well understood. In particular, lepton identification or misidentification probabilities 
are not perfectly reproduced, and this leads in general to an excess of simulated events 
with respect to real data, mainly in the low transverse momentum region. To correct 
for this, rejection factors are computed for each category of identified leptons in the 
simulation, and applied randomly. 

Finally, it is not obvious that lepton tagging calibrated on Z data will produce efficient 
discrimination in the context of Higgs searches at LEP2, especially in 4-jct events. This 
is because in the Higgs search events are reconstructed by constraining the number of jets 
rather than fixing the jet algorithm resolution parameter as was done for the calibration; 
as a result the measured transverse momentum distributions are different. Figure 18 
shows the distributions of the flavour ratios (i.e. the discriminating power) after the 
applied corrections. The data at the Z and at high energy are in fact in approximate 
agreement. 



4.5 Event 6- tagging 

6-hadrons are almost always produced in pairs and the presence of the second 6-hadron 
significantly improves the 6-tagging of an event as a whole. The likehhood ratio method 
provides a simple way for combining the information from the two 6-hadrons. Keeping in 
mind that each flavour is produced independently from all other flavours, one can write 
the equation for the event tagging variable, where the two jets^ are of categories a and 

as: 

It was found, however, that a simpler way of combining the information from two jets 
into a single tagging variable: 

yl} = y.-yp (21) 

works equally well. Here |/q,, yp are given by eqn. (18). The difference between these two 
equations is that (21) neglects the correlated production of the same background flavours 
in an event. The 6-tagging variable computed from (21) was used in the Higgs search 
(see Section 6.5). 

*In events with more than 2 jets, the smallest two values of y were used. 
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Figure 15 compares the performance of event tagging and of jet tagging for Z events. 
As can be seen, very strong suppression of background (down to 10~^ level) can be 
achieved with event tagging. 

4.6 Equalised Tagging 

Physics analyses at the edge of detector capabihties, hke the search for a Higgs boson, 
demand extremely high performance of the 6-tagging. The only way to achieve this 
objective is to expand the set of discriminating variables. But adding a new variable in 
the combining relation (16) becomes more and more difficult with the growth of their 
number due to the increasing influence of correlations among them. 

However, the combined method can be modified to include correlated variables. The 
main idea of the combined method which guarantees optimal tagging for non-corrclatcd 
variables consists in assigning the same value of the tagging variable to different events 
having the same likelihood ratio for background to signal. As described below, the con- 
sistent application of this principle while extending the set of discriminating variables 
gives a desirable improvement of the tagging performance. The simplest way to explain 
this approach is to consider a particular example, obtained from simulation samples and 
presented in figure 19. 

The upper plot in fig. 19 shows the simulated distributions of charged multiplicity N^.^ 
in 6-jets from the process e+e~ — > HZ and in light quark jets from e+e" — > The 
latter process presents the main background for the SM Higgs search, and is suppressed 
mainly by the 6-tagging. The fact that these two distributions are different implies that 
it is useful to include this variable in the tagging. The lower plot in fig. 19 shows the 
ratio RCW^W' /HZ) of the number of fight quark jets from the e'^e~ W^W~ process 
to that of 6-jets from e'*"e~ — > HZ in the simulation, as a function of the tagging variable 

Xj^t = - logio Vjet, (22) 

where yjet is defined by equations (16-17). Two subsamples of events with Nch < 7 
and Nch > 19 are considered separately. As can be clearly seen, for the same value of 
Xjet the ratio R{W^W~ / HZ) in the two subsamples is different. Events with the same 
value Xjet are thus not equivalent; in one subsample they will contain more background 
contamination than in the other. To restore their equivalence, the variable Xj^t should be 
modified in such a way that all events with the same value of X'^^t different subsamples 
will have the same ratio R{W^W~ /HZ). Due to the equivalence principle formulated 
above, such a modification should give better tagging. Technically, equalising of Xjet is 
achieved by a linear transformation: 

X'^ct = A ■ Xjct + B, (23) 

assuming that the dependence of RiyV^W" / H Z) on X^et in each case can be approx- 
imated by an exponential, as shown in fig. 19. The coefficients A and B are different 
for each subsample; their calculation using the parameters of the exponential functions 

is straightforward. 

Including a new independent variable Xnew in the tagging using (16) is equivalent to 
the transformation X'^^t = ^jet - logioyneiy, where ynew = f^^'^i^new)/ f^''^{xnew), i-e. it is 
a particular case of (23). Such a simple transformation cannot be used for Nch because 
of its strong correlation with other discriminating variables, which is reflected in the 
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significantly different slopes of the lines in fig. 19. Instead, the transformation (23) 
works reasonably well. 

A practical application of the equalising method is to the Standard Model Higgs boson 
search. A set of additional discriminating variables is defined for each jet of the event. 
For each new variable, jets are classified in 3 to 5 subsamples. For example, for Nch 
these subsamples are: N^h < 7] 7 < N^h < 12; 12 < N^h < 20; N^h > 20. For 
each subsample the transformation (23) is applied independently and the new tagging 
variable X'j^^ is computed. The parameters of the transformation arc determined from 
the condition that the dependence of RiW^W^ / H Z) on the modified X'-^^ becomes the 
same for all subsamples. The variables are included in the tagging sequentially. For each 
new variable, the X'j^^ obtained at the previous step is used. As before, the global event 
6-tagging variable X^vt is defined as the sum of the two highest X'^^^ values among all jets 
in the event. 

The additional variables included in the 6-tagging using this equalising method refiect 
mainly kinematic properties of 6-quarks. They are: the polar angle of the jet direction; 
the jet energy and invariant mass; the charged multiplicity of the jet; the angle to the 
nearest jet direction; and the number of particles with negative IP. 

Returning to the example of N^h^ fig- 20 shows that equalising over this variable 
improves the suppression of the e^e~ — > W^W~ background. As can be seen from 
fig. 19, the largest difference between subsamples with different Nch is observed at low 
Xjet values, while the background suppression at high Xj^t is almost the same. The main 
improvement from the equalising procedure can thus be expected for the low purity / high 
efficiency tagging, corresponding to low X^et values. Exactly such behaviour is observed in 
fig. 20: including N^h gives almost no improvement for the region of strong background 
suppression. However, equalising the 6-tagging for the complete set of variables given 
above suppresses the e+e~ W^W~ background by an extra factor of more than 2 over 
a wide range of e^e^ HZ efficiency. This additional suppression is important for the 
Higgs boson search since it results in a sizable increase of its detection sensitivity. 

The same equalisation procedure was applied for the e'^e~ — > hA channel when both 
Higgs bosons decay into 66, which is the dominant channel with BR larger than 90% at 
LEP2 energies. A new X^^ was constructed with the condition that all events with the 
same value of X'^^ in different subsamples will have the same ratio R{W^W~ /hA). The 
X^^ variable was used in the search of the hA 4b channel as described in section 6.5. 

4.7 Different Ways of Combining Variables 

As can thus be seen, for combining the separate variables that are relevant for h- 
tagging, three different methods have been used : 

• The IPs of the different tracks are combined by constructing the lifetime probability 
P"*" from the probabilities for the significance values of the various tracks (see Section 
3.5); 

• For the different variables contributing to the 'combined tag', a likehhood ratio 

method is used, as described in Section 4.1; 

• For extra variables, the 'equalisation' method of Section 4.6 is used. 

To some extent these differences are a result of the historical evolution of our h- 
tagging algorithms, but there is some imdcrlying logic to these differences. Thus the 
likelihood ratio method is guaranteed to give the optimal signal/background ratio even for 
correlated variables (assuming of course that the simulation accurately describes the data. 
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including the correlations). However, the method is much simpler when the variables are 
uncorrelated, and this is how the likelihood ratio method was used in 'combined tagging'. 

This could have been used for combining the individual IPs, since the error corre- 
lations between tracks are small. However, it would have been necessary to produce 
signal/background probability ratios separately for each class of track (i.e. for each pat- 
tern of hits in the VD). We preferred to use the lifetime probability, where instead the 
tuning was performed separately for the different track classes. 

Finally, for the Higgs search, extra variables were included to improve the 6-tagging 
performance. Some of these had significant correlations with those already used, so they 
could not simply be added as extra variables in an extended combined tagging approach. 
This led instead to the 'equahsed tagging'. 

5 Modelling and Tuning of Mass-related Parameters 

The agreement between data and simulation is sensitive to the modelling of the physics 
in the event generator, and to the tuning of its parameters. A detailed description of the 
physics model in the main generator used by DELPHI at LEP can be found in [16]. The 
strategy adopted for the parameter tuning and its corresponding results can be found 
in [17]. In this section, two aspects of the modelling are mentioned, which are specially 
relevant in the context of 6-tagging, and which have been investigated recently [18, 19]: 
the modelling of the rate of gluons radiated off 6-quarks relative to light quarks, and 
the probability of secondary 6-quark production through gluon splitting. Both quantities 
affect the description of the dependence of Rb on the jet multiplicity (and on event shape 
variables), and depend critically on the way quark mass parameters are introduced in the 
treatment of the quark fragmentation and subsequent hadronization used in the generator. 
They are important for several of the measurements in Section 6, particularly those 
involving the analysis of multi-jet 6-tagged events, such as the measurement performed 
at LEPl of the running b quark mass at the Mz energy scale (see Section 6.4), and the 
LEP2 Higgs boson search (see Section 6.5). 

5.1 Treatment of Gluon Radiation off 6-quarks 

Discrepancies between simulation and data were observed, which could be attributed, 
entirely or at least partly, to imperfect modelling of mass effects in the generator. As 
an example, the Rh fraction evaluated separately for two and three jet events, using the 
method described in Section 6.1, is illustrated in figs. 21, where JETSET version 7.4 [16] 
was used for the simulation. Similar behaviour is observed in the comparison of two and 
four jet events. 

Because of their higher mass, 6-quarks radiate fewer ghions than lighter flavours. This 
results in fewer multi-jet events in the case of 6-quarks. From kinematic arguments, the 
suppression scales approximately as m^/(s ■ y), where mi,, s and y are the b mass, the 
square of the collision energy and the jet resolution parameter, respectively [20]. It is 
observed explicitly in the value of R'^^, the double ratio of the 3-jet rate for b and light 
quarks,^ used in the measurement of the running 6-quark mass at the Mz energy scale 
(see Section 6.4). Quantitatively the suppression is of order 5%. 

In the original JETSET prescription, used up to version 7.3, mass effects were ignored al- 
together, both in the parton shower evolution describing the fragmentation of the quarks, 

^ilg'' is defined as the ratio of the 6-quark and light quark rates in 3-jet events, divided by the corresponding fiavour 
ratio for events with any number of jets 
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and in the 3-parton matrix element used to correct the first emissions of quarks and anti- 
quarks in the shower. The phase space treatment did include masses, however, and 
induced a large suppression of radiation from the b quark. In version 7.4, and later in 
PYTHIA versions up to 6.130, an intermediate "improvement" was introduced, in that 
matrix element expressions incorporating quark masses were now used in the matching 
procedure. The suppression of the radiation resulting from this intermediate treatment, 
which was in place during much of the LEP period, was however exaggerated by as much 
as a factor of 2, and resulted in the largest discrepancy with the data [18, 19]. Starting 
with PYTHIA version 6.130, and up to version 6.152, mass effects were also introduced in 
the shower evolution through a correction to the expressions of the probabilities of the 
first branchings of each quark. Prom PYTHIA version 6.153, a fully consistent treatment 
is available, including a massive treatment of all branchings in the shower, now taking 
into account in the specification of the matrix element the nature of the couplings of the 
source (vector, axial,...) decaying into quarks, as well as the possibility of unequal quark 
masses (as in the case of W cs). Considerable overall improvement was achieved in 
the description of both 6-tagged 3- and 4-jet rates, thanks to these developments [18, 19]. 

5.2 Treatment of Gluon Splitting to 6-quark Pairs 

Secondary 6-quark pair production from gluon splitting can also result in 6-tagged 
multi-jet events. The corresponding rate is small but is poorly known both theoretically 
and experimentally. This implies an uncertainty in the predictions, particularly of the 
6-tagged 4-jet rate at LEP2 energies. 

Measurements by LEP and SLD collaborations at ^/s = Mz give gyj, = (0.254 ± 
0.051)% [21], where g},i is defined as the fraction of hadronic events containing a gluon 
splitting to a bb pair. This is consistent with the best theoretical estimates, which are 
around 0.2% [18], with relative imcertainties due to unknown sub-leading logarithmic 
corrections, which may be as large as 30%. 

The rate predicted by Monte Carlo generators based on parton shower methods is also 
sensitive to the treatment of sub-leading and kinematic effects in the shower evolution. 
While the original JETSET and PYTHIA prescription resulted in only 0.15%, since version 
6.131 a set of new options has been introduced which bring this rate closer to the measured 
values [18]. Two of these options, which almost exactly double the original rate, have been 
recommended [18] and are used in the latest simulations at LEP2 energies^°. The first 
of these options (MSTJ(44)=3) uses the mass of the virtual gluon involved in a splitting 
to define the scale m^/4 relevant for as, the strong coupling constant, rather than the 
default prescription used for other types of branchings in the shower evolution. The 
second option (MSTJ(42)=3) reduces the conditions on coherence in the emissions, in the 
case of gluon splittings into heavy quark pairs, by introducing a mass correction into the 
angular criterion used to restrict the successive branchings in the shower. 

The impact on the 6-tagged 4-jet rate at LEP2 is best illustrated in the context of the 
Higgs search or of the measurement of Z boson pair production. At LEP2 energies the 
rate of gluons which can split into b quark pairs is less suppressed by kinematics than 
at LEPl. For instance at 189 GeV it is as large as 0.4% using the original JETSET and 
PYTHIA prescription. In a subsamplc of events enriched with 4-jct events, it can reach 
levels near 1% depending on the criterion used on the jet resolution parameter. With 
the new options described above, these values are roughly doubled. The effect of this 
doul)liug on the Higgs soaix li aa ms studied by comparing the numbers of events predicted 

'^"For simulations performed with versions of PYTHIA prior to 6.131, a re-weighting procedure was used to increase the 
g —* bb rate by a factor of about 2. 
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to be selected, when assuming the default value for the rate of gluon splittings into b 
quark pairs or the doubled one. The relative difference between these two numbers is 
about 2.5%; it varies shghtly with the 6-tagging cut but does not exceed 5%. It was 
taken into account in the final evaluation in this channel (see Section 6.5 and references 
therein) . 

6 Physics Applications 

In this section, some analyses involving 6-tagging are described. The aim is not to 
present the physics results, which have already been published, but rather to illustrate 

how 6-tagging works in practice. The extent to which uncertainties in tagging influ- 
ence the final results is also mentioned. Much more detail, including the estimation of 
systematics, can be found in the published papers. 

6.1 Measurement of Ri, at the Z 

One of the most challenging measurements at LEPl is the determination of Rh, the 
branching ratio for Z hadronic decays into 6-quarks. All accurate measurements of Ri, 
use the so-called double tag method. This compares the fractions of events in which 
there is a 6-tag in a single hemisphere with those in which both hemispheres are tagged. 
It allows the extraction from the data of both R^ and the efficiency for tagging a 
hemisphere as coming from a 6-quark. We thus do not have to rely on the simulation for 
the calculation of this important quantity. The small background mis-tag rates, Cc and 
€uds, and the correlation between hemispheres, pb, are taken from the simulation. The 
tracks are separated into two hemispheres by the plane perpendicular to the thrust axis. 
The highest 6-tagging value Xj^t = —logi^yjet (see eqn. (18)) of any jet in the given 
hemisphere is taken as the hemisphere tag. 

The correlation pi, allows for the fact that there are small differences between the 
overall hemisphere 6-tagging efficiency, and the efficiency for tagging a hemisphere if the 
other one has already been tagged. These correlations arise, for example, from the fact 
that at the Z, hadronic events tend to consist of back-to-back jets; if one jet is at large 
positive cos^ where the tagging efficiency is lower (see fig. 3), then the other jet is hkely to 
be at large negative cos^^. again with lower efficiency; thus the efficiencies are correlated. 
The systematic uncertainty on p^ was estimated by comparing data and simulation for 
various kinematic variables that were sensitive to the separate contributions to pi,. 

In particular, it was found that a large correlation arose from the use of a common 
PV for the whole event; if the PV was badly measured as being closer to the SV in one 
hemisphere, then the IP values in that hemisphere would all be systematically reduced 
in magnitude, while those in the opposite hemisphere would be increased. This was over- 
come by determining a separate PV for each hemisphere of the event. This modification 
slightly reduced the flavour discrimination power of the algorithm, and hence increased 
the statistical error, in exchange for a large decrease in the correlation and a smaller 
systematic error. 

Achieving high accuracy for Rh requires the following: 

• the 6-tag must reach very high efficiency to reduce the statistical error: SR^ ~ l/cfe; 

• at the same time the 6-tag must have high purity to reduce the systematic errors 
coming from our knowledge of the background: 5Rb ~ CxRx/^b, where x — q or c; 
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• there must be excellent agreement between data and simulation to reduce the sys- 
tematic errors due to the modelling of the detector resolution, and because there are 
quantities taken from the simulation and not measured directly in the data. 

In the DELPHI Rf, measurement, both the crucial high-purity 6-tag in the "multivari- 
ate" analysis [14] , which finally gave the best precision, and the "combined 6-tag" analysis 
used the combined hemisphere tag described earlier (see Section 4). This required the 
presence of a SV and included the hemisphere lifetime probability Pj, the SV mass Mg, 
the charged energy fraction X^^, and the rapidities of the tracks at the SV. The missing 
transverse momentum P* and the lepton transverse momentum P/ were not used. 

The sources of the systematic error are our knowledge of 6-hadron, charm- and light- 
quark physics (such as lifetimes, decay modes and multiplicities) and our understanding of 
the detector resolution. The first contribution is minimised by measuring the 6-efficiency 
from the data itself and by reducing the charm- and light-quark mis-tag rates to the 
minimum possible so as to give a very pure 6-tag. The second contribution is minimised 
by having good agreement between data and simulation for the detector resolution. 

The high statistical precision of the result is mainly due to the high performance of 
the DELPHI 6-tag: at 98.5% hemisphere 6-purity, the hemisphere 5-efficiency is 29.6%, 
while the mis-tag efficiencies €c — 0.4% and €uds — 0.05%. 

The smallncss of the systematic error comes specifically from the fact that the con- 
tribution from the detector resolution understanding is very small. First, the DELPHI 
vertex detector has three layers of silicon detectors allowing better pattern recognition 
than for a detector with only two layers. The design of the detector is such that the 
intrinsic IP resolutions in both the R(/) and Rz components are good - 27 and 39 iim 
respectively - and consequently also the precision of the primary and secondary vertex 
positions (sec Sections 3.2 and 4.2). Secondly, the detector IP resolution is tuned with 
high accuracy, as described in section 3.6, resulting in a good agreement between data and 
simulation (see fig. 13(b)). As a consequence the error coming from our understanding 
of the detector resolution amounts to only 20% of the total error. 

In summary the performance of the 6-tag and the understanding of the detector resolu- 
tion result in very good stability of the Rf, measurement as a function of the 6-efficiency as 
shown in fig. 22; the highest and lowest efficiencies shown of 44.0% and 21.0% correspond 
to 6-purities of 9L6% and 99.4% respectively. Thus, a stable Rb result was obtained while 
the background contribution varied by more than a factor of 10. The total relative error 
was only 0.4%. 

6.2 R45, the Rate of Events with 4 b-quarks at the Z 

Four b jets are produced predominantly when, in an event with a b pair, a gluon is 
radiated from one of the quarks and itself produces another b pair. This analysis thus 
gives information on the gbb coupling. 

The high purity and efficiency of the tagging method, together with the good agree- 
ment between data and simulation, allowed DELPHI to measure for the first time the 
rate of Z events with 4 6-quarks in the final state [22]. R^b has been measured to be (6.0 
± 1.9 (stat) ± 1.4 (syst)) x 10^^ at a signal efficiency e4b=(3.16± 0.11)%. The analysis 
required 3 jets identified as coming from 6-quarks. This means that the 6-efficiency enters 
to the third power, demanding a very high efficiency of a tight 6-tag. The high purity of 
the tag is required in order to suppress the background from gluon splitting into c-quarks, 
that is a factor of 10 higher than the splitting into 6-quarks. Thus, the measurement re- 
lies on the fe-tagging performance and data/MC agreement in the 6-tagging of the third 
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jet, sorted by decreasing 6-tag value, see fig. 23. The uncertainty coming from tlie 6-tag 
amounts to only 6% of the total systematic error and is determined mainly by the IP 
resolution description. 

This analysis used the 6-tag algorithm of Section 6.1, except that if no SV was found, 
the jet lifetime probability Pj was used by itself. 

6.3 Measurement of the 6-hadron Charged MultipUcity 

The good agreement between data and simulation achieved by the tuning of the track 
resolution allowed DELPHI to measure with very high precision the charged multiplicity 
ub oi weakly decaying 6-hadrons [23]. The basis for the measurement is the determination 
of the number of tracks in a 6-jet which come from the SV rather than from the PV. In 
the hemisphere opposite to the 6-tagged one, the difference 

N+_ ^n-^ -n- (24) 

is computed, in which n"*" and n~ are the numbers of tracks with positive and with 
negative IP respectively. The quantity A^+_ is obtained as a function of the ^-tagging 
purity and the value of is extracted by comparing from data and simulation, 
extrapolated to the limit of 100% 6-purity. 

The result was ub = 4.97 ± 0.03 ± 0.06. The measurement reaches 1.3% precision, due 
to the good understanding of the IP resolution, to the efficient method for determining 
the sign of the IP (see Section 3.3) and to the precision of the VD alignment (see section 
2.2). The tracking efficiency is 99±l%(syst), its uncertainty dominates the systematic 
error on n^. 

This analysis used the same 6-tag algorithm as the previous analysis. 

6.4 Measurement of the Running b-quark Mass at Mz 

The b quark mass determination at the Mz scale has been performed by DELPHI 
by measuring the i?3^ observable, as defined in Section 5.1. Two different jet-finding 
algorithms, Durham [24] and Cambridge [25], were used to reconstruct the jets. Special 
features of this analysis in connection with the flavour tagging performance of DELPHI 
are: 

• b and q initiated events are selected using the same technique, 

• the efficiency versus purity working points are chosen so as to minimise the total 
error on the result (including effects from corrections and biases). 

The directly measured observable, ^ j^ad to be corrected for detector accep- 

tance effects, for kinematic biases introduced by the tagging procedure and for the 
hadronization process in order to get the observable at the parton level, i?3^~^"^ [26]. 
This quantity was then compared with theoretical predictions based on 'Next to Leading 
Order' analytic calculations [20] to evaluate the 6-quark mass at the Mz energy scale. 

In the original analysis [26], q- and 6-quark initiated events were selected by the 
lifetime-signed IPs of charged particles in the event (see Section 3.5). In more recent 
versions of the analysis, the combined tagging technique of Section 4 has also been used. 
The flavour composition of the samples tagged as ^-quark and as 6-quark by the two 
methods are shown in Table 1. 

The magnitude of the corrections applied in each of the two 6-tagging techniques is 
illustrated in fig. 24, where the corresponding ]^^~"^'^°-'^ observables are shown, using 
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Method 


Tagged Sample 


Actually q (%) 


Actually c (%) 


Actually b (%) 


IP 


q 


85.8 


12.6 


1.7 


IP 


b 


5.0 


15.3 


79.7 


Comb 




82.0 


15.5 


2.5 


Comb 


b 


4.3 


10.4 


85.4 



Table 1: Flavour compositions of the samples tagged as g-quark and 6-quark events for 
each tagging method. 

simulated DELPHI data and the DURHAM jet finding algorithm, together with the 
parton-level one, {^^^^""^ obtained with PYTHIA 6.131. For Uc = 0.02, the corrections are 
10% for both techniques, although in opposite directions. 

For the simulation the same corrected result is (obviously) obtained, but for real data 
this is not the case. This difference between the two techniques arises mainly from the 
imperfect modelling of the physics processes affecting the fragmentation and decay of 
heavy and light quarks. Half of this difference is taken as the systematic uncertainty 
associated with this measurement. The error induced by the tagging uncertainties is 
in the range 0.3% to 0.4%, compared with the total error of order 1% on the fiavour 
independence of the strong coupling constant. In terms of the running 6-quark mass, 
these errors correspond to 150 and 500 MeV respectively. 

6.5 Higgs Searches in 4-jet Topology 

One of the main topics at LEP2 energies has been the search for the SM Higgs boson, 
both in the Standard Model (SM) and in the Minimal Super-Symmetric Model (MSSM). 
Here the approach used in the dominant 4-jet channel is outlined. 

Radiative production from a virtual Z e~^e~ ^ Z* HZ is (in principle!) the main 
Higgs process at LEP2 and is referred to as Higgsstrahlung. The mass of the Higgs boson 
is at present unknown, but for a given mass its other properties are determined from the 
SM. The Higgs boson couples to massive fermions and to the and Z bosons. 

The predominant decay mode for the SM Higgs in the mass range of interest for LEP2 
searches is expected to be to pairs of b quarks with a branching ratio ranging from 87% 
to 80% with increasing Higgs mass. The same decay modes are dominant for h and A for 
many choices of values of the parameters in the MSSM, in particular for tan/3 > 1. The 
identification of b jets and rejection of non-6 jets is the most important ingredient in the 
majority of analyses designed to search for the neutral Higgs boson. 

The four-jet final state includes the Higgs boson decaying to 66, and also in principle 
to qq or gg. It is characterised by a large amount of visible energy. As the Higgs boson 
decays mainly to 6 quarks, when the Z also decays to 66 the event will contain 4 6 jets. 
If instead the Z decays to other quark pairs, the topology will again be 4 jets, but with 
only two of them due to 6's. 

The main backgrounds are two-fermion processes e^e~ — > qq{'y), and four fermion 
processes involving VF+VF" and ZZ. Pair production of can result in c jets, but 
only very rarely in 6 jets. The ZZ is an irreducible background if the masses of the 
Higgs boson and of the Z are close, and the Z decays to 66. The cross-sections of the 
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background processes are much higher than that for the Higgs production. At the highest 
LEP2 energies (around 208 GeV) and for uih = 114 GeV/c^, anqq = 0.072 pb while 
(74/ = 19 pb and (72/ = 78 pb. 

In searches for the Higgs bosons at LEP2 in DELPHI, the various differences between 
b jets and light quark jets were accumulated in a single variable Xjet defined for each 
jet, as described in Section 4. Extra variables which help discriminate between the 
signal and background were included in the construction of the equalised tag (see Section 
4.6). Including these extra variables in the tagging algorithm significantly improves the 
rejection of the light quark background. 

The 6-tagging value Xev of the event in the search for the SM Higgs boson is defined 
as the maximum 6-tagging value for any di-jet in the event, computed as the sum of the 
corresponding jet 6-tagging values. In figure 25 the distribution of this equalised 6-tagging 
variable is shown after the common four jet preselection [27] where good agreement 
between data and background simulation can be observed. 

In the top part of figure 26 the performance of combined and equalised methods are 
compared for a SM Higgs boson mass of 114 GeV/c^. As an example using only the 
6-tagging variable, for a signal efficiency of 40% the ^'^(7) contribution is reduced to 
5.6% using combined 6-tagging, wheras it is even more strongly suppressed to 4.3% using 
equalised 6-tagging. The WW component is reduced twice as much using equalised as 
compared with combined 6-tagging. 

In the bottom part of figure 26, the performance of the hA equafised 6-tagging is shown 
in the hA 4b channel and compared with the performance of combined 6-tagging. The 
presence of four 6-jets in the signal makes the analysis of the hA 4b channel different 
from the HZ case, where in most of the cases only two 6-jets are present. The event 
6-tagging value X^^ for the Higgs boson search in the HA — > 4b channel is defined as the 
minimum fo-tagging value for any di-jet in the event and is computed as the sum of the 
corresponding jet 6-tagging values. As can be seen, the application of the hA equalised 
b-tagging improves significantly the performance of the hA — > 4b channel selection. 

As an example using only the X^^ 6-tagging variable, for a signal efficiency of 50% 
the qqij) mistag rate is 0.8% using combined 6-tagging, while it is reduced to 0.5% using 
equalised 6-tagging. The WW efficiency is reduced from 0.1% to 0.06% when changing 
from combined to equalised tagging. The ZZ efficiency is also reduced from 3.7% to 
3.3%. 

In the search for the SM Higgs boson in the four jet channel, the 6-tagging variable was 
combined with another set of discriminant variables [27] using an artificial neural network. 
The final confidence level estimation is calculated using two-dimensional information, 
where one dimension is the neural network output and the other is the reconstructed 
Higgs boson mass. 

6-tagging is also used in the selection of the Higgs di-jet [27] from among the 6 possible 
Higgs di-jet candidates in a 4-jet event. The proportion of correct matchings for the Higgs 
di-jet, estimated in simulated signal events with 114 GeV/c^ mass, is around 53% at pre- 
selection level, increasing to above 70% after the tight selection cut (see table 1 of [27]), 
while keeping a low rate of wrong pairings for ZZ background events. 

7 Conclusion 

The standard approach used by DELPHI for tagging 6-hadrons has been described. 
By using not only the track impact parameters, which are sensitive to the longer lifetimes 
of hadrons containing 6-quarks, but also other kinematic information related to secondary 
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vertices, track rapidities and any leptons, the efficiency/purity has been improved. For 
Z events, a purity of 98.5% for 6-jets was achieved for an efficiency of 30%. Such high 
purity was required for the accurate measurement of R\, at the Z. 

The tagging algorithm can also be applied to complete events, rather than just to 
single hemispheres. For the SM Higgs search at LEP2, the sum of the two largest 6-tag 
variables for jets in the event was used. High efficiency 6-tagging was required in order to 
extract any possible — > 66 decays from the large backgrounds. For a signal efficiency 
of 60%, a rejection factor of 140 for the W~^W~ background was achieved. 

In contrast, the algorithm could be used in an anti-tagging mode, to select jets from 
light-quarks. This was used (together with conventional 6-tagging), to compare the 3-jet 
rates for 6- and for light quarks. This is sensitive to the 6-quark mass. 

For these and for other physics processes, it was crucial to have an efficient, well 
understood procedure for tagging 6 quarks. This resulted in systematic errors being kept 
to a minimum, and enabled many significant measurements to be performed. 
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Figure 1: Schematic cross sections of the Double Sided Vertex Detector in (a) the trans- 
verse view and (b) a three-dimensional view. Dimensions are in cm. 



29 



Outer Layer 
R=106 mm 




Figure 2: Schematic view of the Sihcon Tracker 
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Figure 3: Selection efficiency of Z bb events versus cos Qthmst using the combined event 
6-tagging variable X^v > for the Silicon Tracker and the Double Sided Vertex Detector. 
{Xev is defined in Section 4.1.) The extra coverage provided by the Silicon Tracker in the 
forward direction is clearly visible. 
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Figure 4: The Rep IP uncertainty as a function of p sin 2 (9 (upper plot), the Rz IP uncer- 
tainty as a function of p sin 2^ (middle plot) and the Rz IP uncertainty as a function of 
p for tracks with 6 = [80° : 100°], i.e. perpendicular to the beam direction (lower plot). 
The data are from the Double Sided Vertex Detector at the Z. The curves are parame- 
terisations of a constant intrinsic resolution term and a momentum-dependent multiple 
scattering contribution. Momenta are in Gev/c. 
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R(t) component Rz component 




Pq / ' IPc- ^1 cos e 

IR-RI = (u.V)/sin0 



Figure 5: Definition of R(j) and Rz IP components, m is a unit vector along tlie track 
direction, and ~e is anotlier unit vector in tlie R(j) plane, perpendicular to li. is a 
vector from the origin O to the primary vertex PV. Pq and Pq are the points of closest 
approach in the Rep plane of the track trajectory to O and to PV respectively. The 
diagrams show the projections onto the R(j) and Rz planes. The IP components are dn^ 
and dfiz, while Er^ and Erz are the corresponding components from Pq to the origin. 
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Figure 6: Distributions of the difference between reconstructed and generated directions 
of 6-hadrons for events with a reconstructed secondary vertex, for simulated Z events. 
The B direction is defined as the direction from the primary to the secondary vertex 
(sohd hne) or as the jet direction (dashed hne). Especially in 0, the definition using the 
secondary vertex gives a better description of the B direction. 
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Figure 7: Distributions of positive and negative Rip significances in data Z hadronic 
decays. The excess of tracks with large positive significance is due to long-hved particles. 
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Figure 8: Simulated distributions at the Z of the Rcf) track IP probabilities for different 
quark flavours, for tracks with positive lifetime-sign IPs. 
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Figure 9: Distributions of positive and negative lifetime probabilities {P^ and re- 
spectively) for simulated Z hadronic events. 
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Figure 10: Background suppression in simulated Z hadronic events using lifetime tagging. 
The efficiency is the fraction of 6-jets that are tagged as coming from 6-quarks, while the 
purity is the fraction of tagged jets that are really from 6-quarks. 
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Figure 11: a) and b) The Rcf) and Rz significance distributions for tracks with negative 
IP. The points with errors are real data, the histogram is simulation, c) and d) The 
ratios of these distributions (data divided by simulation) . 
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Figure 12: a) and b) The Rcf) and Rz significance distributions for tracks with negative 
IP after the tuning procedure. The points with errors are real data, the histogram is 
simulation, c) and d) The ratios of these distributions (data divided by simulation). 
The improvement resulting from tuning is clearly visible. 
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Figure 13: The integrated data to simulation ratio of the fraction of selected hemispheres 
as a function of the cut on the 6-tagging variable for a) the lifetime tagging variable 
P^, calculated using all positive lifetime IP tracks in an event (see Section 3.5) and 
b) combined tagging variable Xe^(see Section 4), with tuned (full line) and non-tuned 
(dashed line) track resolution. 
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Figure 14: Distributions of discriminating variables for b- and c- quark jets for simulated 
Z hadronic events, a) Mass of particles in SV. b) Fraction of charged jet energy included 
in SV. c) Transverse momentum at SV. d) Rapidity for each SV track. 
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Figure 15: Background suppression in Z hadronic events using combined 6-tagging. 
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Figure 16: Transverse momentum distribution of identified leptons with respect to the jet 
from which they originate, measured in hadronic Z decays at LEPl. The contributions 
from hght quarks, c-quarks and 6-quarks are added and compared to the data. 
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Figure 17: Distribution of logiQ^P^) for jets with (solid line) and without (dashed line) 
leptons. 
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Figure 18: The upper figure shows the ratio of the numbers of hght-quark and 6-quark 
events as a function of the lepton transverse momentum (dots), as extracted from sim- 
ulated hadronic Z decays. The hne is the fit to that ratio. The lower figure shows the 
same ratio, extracted from a simulated Higgs signal and WW background at LEP2. The 
line is the Z fit. 
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Figure 19: Upper plot: Simulated distributions of the number of tracks in a 6-jet from the 
e+e~ HZ process and in a light quark jet from the e+e~ — > W'^W~ process. Lower 
plot: the ratio R{W^W~ / HZ) of the number of light quark jets from the e^e~ — > W^W~ 
process to that of 6-jets from e^e~ HZ process in the simulation (arbitrary normali- 
sation) as a function of — log^^g Vjeu shown separately for jets with less than 7 or greater 
than 19 tracks. The hnes in each case show the exponential fit of these rates. 



44 



DELPHI 




0.4 0.5 0.6 0.7 0.8 0.9 1 

Efficiency e*e" HZ° 



Figure 20: Mis-tagging efficiency for e~^e~ W^W~ versus e"'"e~ HZ selection 
efficiency, as obtained from standard combined 6-tagging, 6-tagging equalised with respect 
to Nch, and equalised 6-tagging with the complete set of variables (see text for details). 
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Figure 21: Comparison of the measurements of the Rb ratio in 2 and 3 jet events at the 
Z, as a function of the cut in the 6-tagging variable Xj^f The simulation used JETSET 
7.4. 
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Figure 22: Stability of the result as a function of the b-tagging efficiency, for data 
collected in 1994-5. The arrow shows the 6-efficiency chosen for the final result. 
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Figure 23: Distribution of the 6-tagging variables for the first three jets, ordered according 
to their 6-tagging variable, for Z data (dots) and simulation (histogram). The arrows 
show the positions of the cuts used to select the h jets. The figure is from reference [22] . 
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Figure 24: The R'^^ observable from simulated DELPHI data for the two tagging 
techniques and the parton level j^^'P"-"^^ obtained with PYTHIA 6.131. The Durham jet 
finding algorithm is used. 
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Figure 25: Top: distribution of the equalised 6-tagging variable for 4-jet events at 
^/s — 192 — 210 GeV, data (dots) and simulation (solid line). Bottom: ratio of the 
integrated tagging rates in data and simulation as a function of the cut in the equalised 
6-tagging variable. The agreement between data and simulation is satisfactory. 
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Figure 26: Expected SM background mis-tag rates at ^/s = 206.7 GeV, as functions 
of the signal efficiency for a Standard Model Higgs boson of mass 114 GeV/c^ (top) and 
MSSM Higgs boson of mass 90 GeV/c^ and tan/3 = 20 (bottom) when varying the cut 
on the equalised 6-tagging variable (solid lines) and combined 6-tagging variable (dotted 
lines) . 



